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INTRODUCTION 


THE  PRIMARY  responsibility  of  the  air  traffic  controller 
is  to  insure  the  safety  of  the  vehicles  under  his  control. 

He  does  this,  usually,  by  providing  adequate  separation  at 
all  times.  What  is  adequate  separation ?  The  answer  to  this 
question  is  coming  under  increased  scrutiny  as  the  number  of 
vehicles  in  the  air  increases  and  as  efforts  commence  to 
automate  the  controller's  task.  In  terminal  areas,  partic¬ 
ularly,  techniques  for  increasing  airport  capacity,  such  as 
closely-spaced  parallel  runways,  curved  approaches,  speed 
control,  and  V/STOL  operation,  are  limited,  ultimately,  by 
the  closeness  with  which  vehicles  are  allowed  to  approach 
each  other.  Flight  paths  must  be  such  as  to  make  most 
efficient  use  of  the  limited  airspace  available  without 
incurring  an  unacceptable  risk  of  collision.  Very  often 
the  actual  or  proposed  flight  paths  of  two  vehicles  are 
such  that  the  relative  separation  vector  follows  a  compli¬ 
cated  curvilinear  path  with  no  well-defined  point  of  closest 
approach.  A  systematic  design  of  new  procedures,  then, 
requires  a  quantitative  relation  between  any  prescribed  set 
of  flight  paths  and  the  risk  of  collision. 


There  are  at  least  two  ways  to  go  about  assessing  collision 
risk  between  aircraft  in  motion:  worst  case  analysis  and 
statistical  analysis. 

A  typical  worst  case  analysis*  combines  initial  position 
and  velocity  uncertainties,  pilot  control  uncertainties,  and 
possible  changes  in  flight  path  due  to  pilot  or  controller 
intervention  to  determine  the  limits  of  the  airspace  within 
which  the  vehicle  may  be  located  in  the  time  interval  of  interest. 
A  conflict  of  volumes  for  two  aircraft  indicates  that  a  collision 
may  occur;  and  a  warning  is  issued.  This  type  of  analysis  is 
particularly  well  suited  to  a  real-time  collision  warning  system. 
Since  it  makes  no  assumptions  about  pilot  or  controller  inten¬ 
tions,  it  can  handle  VFR  as  well  as  IFR  traffic.  But  because 
great  latitude  must  be  allowed  for  the  intentions  of  VFR  air¬ 
craft,  the  volumes  calculated  expand  rapidly  in  time.  This 
limits  the  period  of  applicability  to  the  order  of  seconds  of 
time.  Therefore  it  is  usually  impractical  to  use  worst-case 
analysis  to  design  most  terminal  or  enroute  traffic  patterns;  its 
major  utility  is  in  detecting  and  resolving  potential  conflicts 
in  real  time. 

The  statistical  analysis  of  collision  risk  assumes  that  the 

pilot  will  proceed  approximately  along  the  path  for  which  he  has 

« 

obtained  clearance  from  the  controller,  and  that  jbhe  controller 
will  not  direct  sudden,  unsafe  maneuvers.  For  these  reasons, 
the  method  applies  well  to  IFR  traffic  over  a  period  of  several 
minutes  or  even  hours.  The  method  has  two  aspects:  (a)  the 
collection  and  analysis  of  data  to  determine  the  statistics  of 
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vehicle  position  and  velocity  deviations  from  the  prescribed 

path,  and  (b)  analysis  of  the  distributions  thus  obtained  to 

give  the  probability  of  collision.  The  first  aspect  requires 

combining  equipment  performance  estimates  with  the  results  of 

2 

pilot-vehicle  simulations  .  Very  often  the  only  information 
obtained  about  the  distributions  are  their  variances.  As  a 
result,  the  second  aspect,  which  this  paper  treats,  usually 
employs  probability  distributions  based  on  analytical  con¬ 
venience  rather  than  fact.  Nevertheless,  a  series  of  analyses 
3, 4, 5, 6,  based  Qn  assumed  distributions,  have  developed  risk 
estimates  for  parallel  airways  and,  recently,  for  parallel 

7 

runways  .  A  large  number  of  cases,  however,  cannot  be  handled 
with  present  theory.  Present  theory  deals  directly  with  air¬ 
craft  flying  parallel,  straight-line  paths,  and  cases  reducible 
to  that.  It  cannot  handle  curved  relative  paths,  ascending  or 
descending  airways,  or  crossing  airways,  although  such  cases 
are  very  common  in  the  critical  terminal  areas  in  which  most 
of  the  present  congestion  occurs.  It  is  the  purpose  of  this 
paper  to  provide  the  theoretical  basis  for  calculating  collision 
probabilities  from  the  statistics  for  any  two  given  paths  in 
space.  It  is  hoped  that,  as  data  are  gathered  for  better 
definition  of  the  distributions,  the  formulas  here  presented 
will  assist  in  the  development  of  improved  terminal  area  traffic 
procedures  and  rules. 


-3- 


STATEMENT  8F  PROBLEM  AND  OF  ASSUMPTIONS 

Assume  that  two  aircraft  start  on  prescribed  flight  paths 
at  time  tg.  Due  to  navigation  and  control  uncertainties  the 
actual  paths  deviate  from  the  prescribed  paths.  The  question 
is:  what  is  the  probability  of  collision  from  time  t q  to  some 
later  time,  t^? 

The  following  assumptions  are  made: 

(1)  All  possible  flight  paths  are  random  samples  from 
ensembles  of  known  statistics;  the  means  of  the 
ensembles  are  the  prescribed  paths. 

(2)  Each  vehicle  may  be  represented  as  occupying  a 
bounded,  closed  region  of  space. 

(3)  Vehicle  position  is  a  differentiable  function  of  time. 

(4)  The  mean  separation  is  always  large  compared  to  the 
diameter  of  the  regions  occupied  by  the  vehicles. 

The  first  assumption  is  the  basis  of  the  statistical- 
probabilistic  analysis  ~f  aircraft  collision  hazards;  its 
advantages  and  disadvantages  already  have  been  discussed. 

The  second  assumption  is  needed  to  define  a  collision 
mathematically.  But  it  is  more  than  a  mathematical  nicety,  for 
by  appropriate  selection  of  the  vehicle  volumes  one  may  take 
account  of  wake  turbulence,  or  calculate  near-misses  or  air¬ 
space  conflicts  instead  of  actual  collisions. 

The  third  assumption  is  justified  by  the  finite  accelera¬ 


tion  of  aircraft. 
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ANALYSIS 

i 

Probability  of  Collision 

■  i  !  : 

j  » 

It  is  desired  to  find ‘the  probability  of  at  least  one 

«  i  •  / 

collision  between  the  two  vehicles  in  the  interval  t^  -  tg. 

i 

(More  than  one  collision,  although  practically  impossible,  is 

i 

not  explicitly  excluded  by  the  assumptions.  To  do  so  would 

t 

make  the  use  of  probability  distributions  extremely  complicated.) 

i  ■  ’ 

Let  'N  be  the  number  of  collisions  between  the  two  vehicles  in 
,  » 

I 

.  time  t^  -  tg.  The  desired  probability  is  P[N  -  1] : 

P[N  >  1]  =  P [N  =  1]  +  P [N  =  2]  +  P[N  =  3]  +  •••  . 

'  '  : 

The  terms  P[(N  =  l)Pi(N  =  2)],  etc.',  do  not  enter  because  N 
cannot  have  two 'values  at  once. 

i  *  ; 

Instead  of  calculating.  P [N  £  1]  directly,  it  is  easier 
to  approximate  it  by  N,  the  average  of  N  taken  over  a  large 
number  of  trials’  which  are  ;Ldentical  except  for  statistical 
fluctuations  in  the  two  paths.  (The  tilde  (~)  above  a  quantity 
, indicates  the  ensemble  average  of  the  quantity?  a  vector  is 

I  I 

indicated  by  a  bar  —  above  the  corresponding  scalar  magnitude. 

i  ! 

The  terms  mean,  average expected  are  used  interchangeably  in 

this  paper  to 'refer  to  the  expectation  value  of  a  random 

•># 

variable.)  By  definition,  Ni  is  the  mean  of  the  probability 

i  i 

..distribution  P[N  =  i] .  That  is,  N  =  0»P[N  =  0]  +  1»P[N  =  1] 

»  * 

+  2»P[N  =  2]  +  3*PpN  =  3]  +  .  The  approximation  to  be 

,  i 

used  is  J? [N  £  1]  ~  N.  The  error  in  the  approximation  is,  with 
P[i]  written  for  P[N  =  i] ,  i  =  1,2,3, , 

!  I  i 
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N  -  P[N  >  1] 


=  ^2  i  P[i]  -  P[i]  =  S  (i  “  1)  P[il 
i=l  i=l  i=l 


=  ^  i  P[i  +  1]  ±2  i  P[i]a 


where 


a  >  P[i  +  1]/P[i] 


N  -  P[N  >  1]  <  Vi  pflja1 


<  ?  i  P[N>l]a^,  since  P[l] 


-  -<—J 

1 


=  P  [N=1J  <  p  [N>1] 


The  percent  error  in  the  approximation  is 

00 

(N  -  P[N  >  1])/P[N  >  1]  <  ^  ia1  =  a  /  (I  -  a)2, 

i=l 

assuming,  a  <  1. 

The  percent  error  is  small  provided  only  that  a  is  small. 
The  ratio  a  is  just  an  upper  bound  on  the  probability  of  N+l 
collisions,  divided  by  the  probability  of  N  collisions.  The 
condition  on  a  is  assured  by  qualifying  assumption  (1)  as 
follows:  The  ensemble  statistics  are  such  that  the  probability 
of  N  collisions  between  the  two  vehicles  in  time  t^  -  t^  is 
large  compared  to  the  probability  of  N+l  collisions,  for  N>0. 

The  problem,  then  is  to  calculate  N,  the  average  number 
of  collisions.  In  order  to  do  this,  it  is  necessary  to  make 
precise  the  notion  of  collision. 


Description  of  "Collision  Surface" 


By  assumption  (2) ,  each  of  the  two  vehicles  may  be  repre¬ 
sented  as  occupying  a  bounded,  closed  region  of  space,  say 
V1  and  V2*  The  center  of  volume  of  is  defined  as  usual  to 
be  rx,  ( 

*/v  idvi 7  j'yj  dv1)  and  so  also  for  and  r^. 

The  position  vectors  and  pj  must  have  a  common  origin. 

The  collision  surface,  S  ,  is  now  defined  as  the  locus  of 

c 

obtained  by  translating  V ^  such  that  V ^  and  have  one  or 
more  common  boundary  points  but  no  common  points  that  are  not 
boundary  points.  In  other  words,  the  collision  surface  is 
the  center  of  the  second  vehicle  when  translated  without  rota¬ 


tion  to  touch  but  not  penetrate  the  first.  A  two  dimensional 
picture  of  a  simple  collision  surface  is  given  in  Fig.  1. 

Note  that  the  orientations  of  Vj_  and  V2,  and  therefore  the 
shape  of  Sc,  may  vary  in  time.  Also,  Sc  is  the  boundary  of  a 
region  Vc,  which  contains  r^  Finally,  it  should  be  noted  that 
the  collision  surface  obtained  by  translating  to  touch  V2 
is  the  reflection  of  Sc  about  the  point  The  subsequent 

analysis  and  the  collision  probabilities  obtained  from  it  are 
insensitive  to  which  collision  surface  is  used. 


A  collision  now  may  be  said  to  occur  whenevar  the  rela¬ 


tive  position  vector  r  enters  Sc.  However,  this  is 

not  completely  satisfactory  if  S  has  concavities,  as  shown  in 

c 

Fig.  2.  For  if  the  surface  has  concavities,  the  trajectory 


may  enter  Sc  more  than  once  in  what  should  be  counted  as  a 
single  collision.  There  are  at  least  three  ways  to  avoid  the 


problem:  (a)  use  a  convex  surface  in  place  of  one  with  conca¬ 
vities,  (b)  use  a  surface  with  concavities,  but  allow  for  the 
over-estimate  in  collision  probability  that  will  thus  result, 
(c)  use  the  collision  cross  section  vector,  to  be  defined 
subsequently.  The  first  of  these  alternatives  is  preferred, 
since  it  greatly  simplifies  the  computation.  But  only  the 
third  is  rigorously  correct. 

Average  Number  of  Collisions  in  -  tg 

It  is  now  possible  to  calculate  N,  the  average  number  of 
collisions  in  time  t^  -  t^.  The  approach  to  be  taken,  follow- 

g 

ing  Rice  ,  is  to  construct  an  expression  for  dN/dt,  calculate 
its  statistical  mean,  dN/dt,  and  then  integrate  in  time  to  get 
N.  The  mathematical  legitimacy  of  the  approach  need  not  be 
discussed  here. 

An  expression  for  dN/dt  may  be  built  upon  the  previous 
description  of  collision',  a  collision  occurs  whenever  the 
relative  position  vector  r  enters  the  collision  surface  Sc. 

At  such  time,  N  increases  abruptly  from  zero  to  unity.  Let 
<p  (r)  be  a  step  function  in  ?  space  that  is  unity  for  ?  in  V 

C 

and  zero  for  ?  outside  V„. 

c 

Then  N  may  be  defined: 

N (t)  =  u(v*n)  d<J>  (?)/dt 

where  v  is  the  relative  velocity  (=dr/dt) ,  ii  is  the  inward 
normal  to  S  at  point  of  entry  and  u (x)  is  zero  for  x  <  0  and 

C 


-10- 


I 


unity  for  x  >  0.  :The  purpose  of  u(v*n)  in  this  expression  is 

■  i  |  * 

to  exclude  exits',  for  which  v*n  <  0,  but  to  “include  entries 

■  ’  ,  i  '  1 

{v*n  >  0)  and  mere  touches  (v*n  =Q).  ! 

The  time  derivative  <j>  (r)  may  be  written  as  d<f>/dr,  which  is 
1  1  !  ! 
a  delta  function,  times  dr/dt  provided  the  position,  r,  is  a  diff 

.  .1  '!  : 

erentiable  function  of  time.  Since  this  is  no,  by  assumption 

*  1  !  V 

(3),  it  follows  that  i 

!  ■  1  '  !  i 

N  =  u(v*n)i  (d<f>/dr)  (dr/dt)  =  u (v»n)v*V<j>. 

i 

I  +*  t 

To  obtain  the  ensemble  average  N,  it  is  necessary  t.6  con¬ 
sider,  instead  of  the  single  relative  trajectory,;  an  ensemble  ■ 
of  trajectories.  At  any  instant,. the  probability  density  > 

:  i 

distribution  W(r,v,t)  of  the  relative  position  and  velocity  is 

■  ’  .  .  , 
known,  by  assumption.  Therefore,  the  ensemble  average  may  bp 

>  •  1  ii 

found  as 

.N  =  J' d(v)  ■J' d  (r)i  W(r,v,t)  u(v»ii)  v*^<|>.  i 

Here  d(v)  and  d(r)  are  elements  of  velocity  space  and  position 
space.  The  integrations  are1  understood  to  cover  the  spaces 

i  -  1  - 

completely.  1  By  utilizing  the  identity  V*  (<J>x)  =  x*Vy  +  <j>y*x, 
where  if  is  a  scalar  and  5  a  vector,  say  x  =!v  u(v*n)  W(r,v,t), 

'  ’  j 

the  above  is  written  ,  ,  1 

(  » 

N  =  /d(v)  J" d  (r)  x*V‘  <f)  !  .  > 

i  .  *  * : 

=  J d(v)  j -d{r)  (4>x)  r  <J>  V*xJ 

.  :  ' 

)  1 

=  /d(,)[  J  d  (r)  V  -  (<f>x)  -  f d(y:)  <f> 
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i  '  .  (  , 

By' the  divergence  theorem  the  first  integral  in  the  brades 

which  is  oyer  the  volume  of  ?  space,  is  converted  to  a  surface 

:  '  1 

integral  over  a  sphere  SR  of  increasingly  large  radius  R,  cen¬ 


tered  a,t  r  =  0: 


/ 


d (r)  V* (4>x) 


=  lim ;  /c 
!  R+oo  J  b 


g  (ds  ♦  (<Jix) . 


As  R-*-°°'  the  surfa'ce  SR  eventually  lies  wholly  outside  of  Sc. 
Since  <j>  =  0  outs,ide  the  limiting  value  is  zero. 

The  pecopd  volume  integral  in  the  brackets  is  zero  in  the 

*  ,  i 

space  outside  S  ,  because  4>  is  zero  there.  It  has  a  non-zero 

i  c  !  ; 

value  only  within  t'he  collision  volume  V  ,,  .where  <p  =  1. 

I  i  ■  c 

Therefore  i  .  '  1  ! » 


N  =  -  d (r)  (1)  V*x 

% 

=  /dtf)  Js  ds'*x, 


where  ds  =  (ds)n 


•  i 


f<Hv)  L  (ds"v)  u(:v*n)  W(rc,v,t) 


N  =‘ 


,/»(v 


(v)  W(0,v,t)v* 

i 

1  t 

N  dt  . 


ri 


I  ds  u(v*n) 


■  i 


T,wo  notes  regarding  the  above'  are  in  order. 


i  <  .'<> 


v„  i 


i 


I 
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First,  rc,  the  relative  position  vector  drawn  to  the 

collision  surface  S  has  been  replaced  by  6,  the  null  vector 

c 

in  r  space.  The  reason  is  that  the  distribution  W  in  r  will 
have  a  mean  large  compared  to  r  in  all  practical  cases,  as 

w 

stated  in  assumption  (4).  The  difference  in  the  density  W  for 

r-  =  r  and  r  =  0,  therefore,  is  small.  As  a  result,  W  may  be 
c 

taken  outside  the  surface  integral. 

Second,  the  surface  integral  ds  u(v*n)  may  be  inter¬ 


preted  as  a 


'ision  cross-section  vector  with  respect  to  v, 


and  written  AS„(v/.  If  has  no  concavities,  the  vector  AS 
c  c  c 

has  the  property  that  its  component  in  the  v  direction  is  the 

area  of  the  orthogonal  projection  of  S  onto  a  plane  perpendi- 

cular  to  v.  If  S_  does  have  concavities,  then  this  property, 

c 

rather  than  the  surface  integral,  serves  to  define  ASc.  in 


either  case  one  has 


" =  fc  [/ 


d(v)  W(0,v,t)  v»AS  (v) 

c 


The  general  formula  thus  derived  is  a  three-dimensional 

form  of  the  well-known  expression  for  the  mean  number  of  times 

a  random  function  y(t)  crosses  a  given  level  yQ  in  a  time 

interval  t^  -  tQ.  The  one-dimensional  formula,  originally 

8 

derived  by  S.O.  Rice  ,  is 


/,  00 

f  dt  /  W1(yQ,y,t)  |y|  dy 
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v/here  w  is  the  joint  probability  density  of  y  and  y.  The 
three-dimensional  formula  gives  the  average  number  of  times  a 
random  vector  r(t)  enters  a  given  closed  surface  S  in  time 

vC 

t.^  -  tQ.  The  major  difference  between  the  two  formulas  is 
that  the  one-dimensional  case  requires  specification  of  a 
level,  Yq,  while  the  three-dimensional  case  requires  specifi¬ 
cation  of  an  entire  surface,  S  .  The  term  |y|  of  the  former 

C 

case  becomes  the  more  complicated  integral  v*Ss"  : 


?-ASc  =  V-  L  u{v*n)  n  ds 


Special  Case:  Large  Relative  Velocity 

The  simplest  case,  mathematically,  occurs  when  the  mean 
relative  velocity  is  large  compared  to  the  spread  in  the  dis¬ 
tribution  of  velocity-  Fortunately  this  case  covers  many 
practical  situations,  such  as  curved  landing  approaches, 
crossing  airways,  climbing  and  holding  patterns  and  over¬ 
flights  of  landings  or  takeoffs.  The  situation  may  be  approx¬ 
imated  by  assuming  the  density  W(0,v,t)  to  be  a  delta  function 
in  velocity  space  at  the  mean  velocity  v.  In  that  case, 


/• 


d  (v)  v  W(0,v,t)  =  v  (t )  Wr(0,t) 


where  W  (r,t)  is  the  density  distribution  in  r.  As  a  result, 

<V 

N  becomes 
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tl 

N  =  f  dt  v  (t)  •  AS  W  <0,t) 

Jtn  c  r 


t, 

- 


dt  (dr/dt)  •  AS  W(0,t) 

C  i 


r(t, ) 

£ 

'r!t0) 


dr  •  ASC  W£(r,t) , 


where  W'  is  W  with  mean  translated  to  0. 
r  r 

In  this  formula,  N  may  be  visualized  as  the  volume  swept 
out  through  the  position  distribution  Wr  by  the  collision  sur- 
face  Sci,  in  going  from  r<tQ)  to  r(t.).  Such  a  picture  is  use¬ 
ful  if  the  relative  path  is  irregular  or  if  the  density  distri 
bution  is  time -dependent.  In  such  cases  it  is  practical  to 
program  the  integration  of  (2)  on  a  digital  computer  using 
relatively  large  step  sizes  and  simple  integration  algorithms. 

*v 

An  order  of  magnitude  estimate  of  N  is  usually  adequate. 

If  the  relative  path  is  a  straight  line  and  if  the  den¬ 
sity  distribution  is  time-independent  and  analytically  conven- 

♦ 

ient,  further  analysis  of  equation  (2)  can  be  carried  out. 

The  following  analysis  assumes  straight  line  paths  and  time- 
invariant  Gaussian  statistics,  a  case  treated  ad  hoo  in  the 
literature1'3'4. 

Let  the  vehicles  be  represented  by  two  right  circular 
cylinders  with  parallel  bases  of  diameters  d^  and  d2  and 
heights  h^  and  h2»  Let  the  mean  paths  r^(t)  and  ?2 (t)  be 
horizontal  straight  lines  going  to  infinity  in  both  directions 
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finally,  let  the  actual  positions  and  r2  be  normally  and 
independently  distributed  with  stationary  statistics  about  the 


mean  positions: 


"rl'V  *  W3'2  lKrl< 


-1/2 


Wr2(r2}  =  (27r) 


-3/2 


-1/2 


exp  -  |  ' 

exp  “  j  ^2"^2^/r2  * 


Here  subscript  T  denotes  transpose  and  K~J  and  K~2  are  the 
inverses  of  the  covariance  matrices  Krl  and  Kr2.  These  co- 
variances  are  assumed  to  be  known  in  the  (x^y^z^)  and  (x2y2z2) 
coordinate  systems  illustrated  in  Fig.  3.  In  the  figure, 

A*  <V 

the  position  errors  (r1  -  r^)  and  (?2  -  r2)  have  components 
x ^  and  x2  in  the  direction  of  travel,  and  z2  along  the 
vertical  and  y^  and  y2  normal  to  the  other  two  components.  In 
these  coordinates,  assuming  the  position  errors  are  statisti¬ 
cally  independent  in  the  three  directions,  the  covariances  are 


Krl  -  0 


To  calculate  N  in  this  case,  it  is  necessary  first  to  com¬ 
pute  the  distribution  of  r2  -  ^  and  then  to  compute  the  por¬ 
tion  of  that  distribution  swept  out  by  a  right  vertical  circu¬ 
lar  cylinder,  height  h.^  +  h2  and  diameter  d-j^  +  d2,  moving  on 

£ t  £ 

the  straight  line  r2  -  r^. 
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The  distribution  of  ?2  -  5^  is  normal#  being  the  sum  of 
normal  random  variables.  If  the  position  errors  of  the  two 
vehicles  are  independent,  as  assumed  for  this  example,  then 


Kr  “  Kr2  +  Krl  = 


2  2  2 

a  a  a 

x  xy  xz 

2  2  2 

a  a  a 

yx  y  yz 

2  2  2 

0  0  0 

zx  zy  z 


The  oz  terms  here  are  determined  by  writing  Krl  and  Kr2  in  the 
same  (x  y  z)  coordinate  system,  which  is  yet  to  be  selected, 
and  adding.  The  distribution  Wr(0)  is 

Wr(0)  =  (2*)~3/2  lKrf1/2  exp  -  |  [(-r)  K"1  (4)J 

s  W^.  (r) . 

The  portion  of  this  distribution  swept  out  by  the  collision 
surface  is  given  by  formula  (2).  Choosing  x  to  be  along  the 
relative  velocity  vector,  y  horizontal  and  perpendicular  to  x, 
and  z  upward,  as  shown  in  Fig.  4,  formula  (2)  leads  to: 

N  =  (h-j+h2)  (d1+d2)/(2Tt  0z0y)  exp  -  1/2  ^(y2/oy2)  +  (z2/az2)j 


where 


jx2  =  a12cos2a  +  a22cos20  +  alysin2a  +  o2ysin2e 
Ty2  =  alycos2a  +  a22cos23  +  o2xsin2a  +  022sin2$ 


2  _ 


xy 


|<°ly  '  °lx)sin  2a  +  ?(c2y  '  °2l)sin  26 
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a 


I 


$'  i 


•  i  ; 


>2  _  2  .  '  2 
°Z  =  °lZ  +  a2z 


a  2  =  a  2  «  0 
xz  yz 


u  =’  a  /  a  a 
p  xy  '  x  y 


For  the  special  case  of  Vehicles  passing  on  the  same  level  in 

\ 

\  1 

opposite  directidns  with  mean  path  separation  S^, 


y  =  Sy' 


Z  W  0 


a  =  it 


3  =  0 


2  _  „  2  ;  2. 
cr  =  0  _  +  a,  , 

y,  2y  ly  ' 


_  2  '2  .  „  2- 
az  =  a2z  +'alz  ' 


5=  ((Vhj)[artl/I.)/1(l.1^;)W  t'’ly«2y)1/2) 

r  i  [<Sy'/<'-  2y+aly>] 


.  eXP 

-  I  J  ~J  “J  J 

i  u  2  2 

•If,  further,  if  is  .assumed  that  h^  =  h2,  d1  =  d2,  o ly  =  c2y, 
and  alz  =  then  the  value  of  N  ‘coincides  exactly  with  that 

obtained  by  Riich3  and  Taylor.4  : 

i  i 

Speiif  Can:  Nan-Zara  Rafativa  Valacity 

If  the  mean  relative  velocity  is  not  large  compared  to 
the  higher  moments,  the  delta-function  approximation  for  the 

•  ,  i 

, velocity  distribution  cannot  be  made.  This  situation  occurs 

.  i  . 

when  there  is  a  small  (but  non-zero)  relative  velocity,  as,  for 
example,  when  two 'aircraft  are  controlled  to  be  in  train  pre- 
•paratory  to  final  approach  ,  or  to  be  following  parallel  en- 
route  airways  but  with  slightly  different  speeds. 

J.t  will  be  assumed  for  rthe  present  case  only  that  the 

mean  relative  velocity  ajid  the  higher  moments  are  all  different 

»  \  ‘ 

from  zero  and  that  the  position  errors  are  independent  of 
velocity  errors.  Then  the  joint  distribution  for  each  vehicle 
may  be  factored;  ' 

l  ; 
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i 


9 


Vri.Vt)  s  wlr(rlft)  Wlv^l'tJ 
w2(?2,v2,t)  =  w2r(J2,t)  w2v(J2,t;  . 

Since  relative  position  r  and  relative  velocity  v  are  sums  of 
random  variables,  they  are  independent  if  (r1+r2)  is  independ¬ 
ent  of  (v^v^.  If  so,  W(?,v,t)  also  may  be  factored  so  that 
formula  (1)  is 


N  -y4  dt  Jd(v)  Wv(v,t)  v-Ssc  Wr(0,t)  .  (3) 

If  the  collision  cross-section  and  the  velocity  distribu¬ 
tion  are  independent  of  time,  the  portion  in  square  brackets 
may  be  taken  outside  the  time  integral.  Nevertheless,  a  com¬ 
puter  solution  is  usually  necessary.  One  case  that  can  be 
handled  analytically  is  that  of  one  aircraft  slowly  overtaking 
and  passing  another  on  a  parallel  airway,  with  normal  statistics. 

9 

The  analysis  is  tedious  and  hence  only  the  result  will  be  given 


here. 

If  the  statistics  are  normal  and  if  the  vehicles  are  repre¬ 
sented  by  vertical  circular  cylinders  of  heights  h^  and  h2, 
diameters  d1  and  d2,  (see  Figure  3)  moving  on  horizontal 
straight  paths,  then  the  collision  probability  is  approximately 


,  /T2 — 72 

(d1+d2)2A^.(d1+d2)  (hj+h^/al^  +A*  f(v,K)/(4 /Sr  o  o  | v | ) 


•  exp  -  j 

m 


(y  2/ohHz2/& 
y  z 


21- 


m  i 


v'  =  v2  sin  (8  +  8)  -  v.^  sin  (a  +  6)  =  v  sin  6 


,  ,2  _  1  M  2  ,  .  2  ,  ,  2,,  2. 

A  x  =  I  <Alx  +  A2x  +  Aly  +  A2y> 


"  H 


<A?y‘AL>2  +  <Aiy-A2x>2 


+  2  <Aly-AL>  (A2y-A2x>  cos  2  <* 
A'y  5  k  (Alx  +  A2x  +  Aly  +  A2y> 


-*0 


1/2 


l(< 


+  -<<-Alx>2  +  <A2y-AL>2 


+  2  (^y-^x)  (^y-AL>  cos  2  <a 


-B)) 


1/2 


ax'ay,0z'°xy,axz'ayz:  as  the  case  of  large  relative  velocity 


X » 2  =  ^2  ,  i  2 

A  z  ~  Alz  +  n2z 


tan  26  =  [(Xiy-A^x)  sin  2a  +  sin  2f3]/[(A 

2  \2 
l2y"A2x 


2 

iy 


-aJx>  cos  2a  +  (*2y-A2x}  cos  2e] 
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S&X^TfflHMar 


2  2 

In  the  above,  vehicle  1  has  position  variances  alx,  aly, 

2  2  2  2 
°lz  velocity  variances  Alx,  A^,  A^z  which  are  along  its 

track  (x) ,  across  its  track  (y) ,  and  along  the  vertical  (z) ;  its 

mean  velocity  v^  makes  an  angle  a  with  the  mean  relative  velocity 

~  2  2  2  2 
vector  v.  (See  Fig.  4)  Analogous  terms  cr2x,  a2y,  a2z,  A^, 

2  2- 

A2y,  A2z,  v2,  3  hold  for  vehicle  2.  The  mean  vertical  path 
separation  is  z,  the  closest  horizontal  approach  of  the  mean 
paths  is  y.  The  six  position  and  velocity  errors  of  a  vehicle 
have  been  assumed  uncorrelated  with  each  other  and  with  those 

of  the  other  vehicle. 

Special  Cate:  Zero  Aelafive  Velocity 

When  two  aircraft  using  the  same  airway,  or  in  train  to 
the  same  landing  approach  are  controlled  so  that  (on  the  average) 
their  separation  is  constant,  then  the  mean  relative  velocity 
is  zero.  This  case  must  be  handled  separately,  because  if  in 

<v 

the  preceding  result,  the  mean  relative  velocity  v  is  allowed 
to  approach  zero,  the  value  of  N  increases  without  bound. 
Physically,  such  a  result  is  expected,  since  N  is  constant  in 

CO 

time,  and  /  dt  N  must  be  unbounded.  The  infinite  limits  of 

I  ‘/-CO 

time  should  be  replaced  by  the  more  realistic  limits  tg  and  t^. 

If  the  statistics  are  stationary  with  velocity  independent  of 

<v 

position,  the  formula  for  N  is 
t. 


N 


-/>  / 


d (v)  W(0,v)v*£s' 


“  (tl  ”  V 


/ 


d(v)  W  (v)v*AS_ 

V  C 


Wr(0) 


(4) 
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To  illustrate  an  analytic  solution  in  the  case  of  zero 

relative  velocity,  the  same  geometry  will  be  assumed  as  for 

the  two  preceding  cases,  i.e.,  two  right  cylinders  moving  on 

the  mean  straight  horizontal  paths  shown  in  Fig.  3.  (The 

statistics  are  assumed  to  be  gaussian  and  stationary,  with  no 

9 

correlation  between  vehicles.)  The  result  obtained  is: 


Sptsiil  Cut:  Sphtriea!  Ctllitiw  Surfact 


When  the  collision  surface  Sc  is  spherical  the  collision 

«•  A  2  A 

cross  section  vector  is  vira  where  v  is  a  unit  vector  in 

c 

the  direction  of  v  and  a(t)  is  the  radius  of  the  sphere  at 
time  t.  Then  the  expression  for  N  is 


/>  /• 


d(v)  W (o, v, t)  v»v  it  a 1 


wr  (o,t)  v(t)  TT  a2  (t) 


(5) 
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where  v(t)  is  a  scalar, -equal  to  the  average  magnitude  of  dr/dt. 

,  i 

This  quantity,  v(t),  is  just  'the  mean  relative  speed;  it  is 
never  negative  and  is ■ zero  only  when  the  velocity  distribution 
has  zero  riean  and.no  spread  (a  delta  function  at: the. origin) .  , 

*  i  ,  • 

Unlike  (2) ,  formula  (5)  is  not  restricted  to  large  relative 

'  !  j 

velocity.  Its  only  limitation  is(that  it  applies  to  spherical 
collision  surfaces.1  , 

In  many' cases. an  actual  collision  surface  pay  be  replaced 

.  .  i 

by  a  spherical  one  of  equal  area  without  seriously  impairing 

.  ■  >  •  ! 

the  usefulness  of  the  result.  This  occurs, i  for  'example,  in 

: 

determining  the  probability  of  near  pisses,  or  ^.n  obtaining 

an  order-of -magnitude  answer.  If  such  an  approximation  may  1 

f  . 

be  made  then  the  collision  probability  is  given  by  (5)  with 

•  ^  *  • 

A (t) ,  one  fourth  the  area  of  the  actual  collision  surface 

2 

at  time  t,  in  place  of  xa  (t) .  1  1 


i  i 
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APPLICATION  TO  PARALLEL  RUNWAYS 

i  1 

•The  use  of  closely  spaced,  independently  operated, 

i  !  " 

>  >  .  .1:  .  i  . 

parallel  runways  is  one  means  to  increase  airport  capacity 

.  ,  .  1 

that  was  recommended  for  further  investigation  by  the  Air 

1  10  .. 

Traffic  Advisory  Committee  Report.  The  scheme  involves 

i  * 

multiple  curved  precision  approaches  to  two  or  more  parallel 

:  ,  I 

'runways  spaced  ;at  less  than  the  present  500d-ft  minimum.  No 
time  synchronism  is  assumed  between  aircraft  ion  different 
tracks.  The  safety  of  such  approaches  idepends  on  at  least 

n  ! 

four  factors:  the  probability  of  missed  approaches,  , 

i  • 

emergency  (e.g.1.,  engine  out)  procedures,  accuracy  of  guidance 

■  1  i  i  .  i 

of  each  aircraft  relative  to  its  runway,  ’accuracy  of  guidance 

,i  i 

of  each  aircraft  relative  to  thesother.  The  last  of  these 

■  i  ' 

factors  &as  analyzed  in  detail  by  the  method  developed  in 
this  paper,  and  that  analysis  will  be  described  next. 

i  i  (  1 

A  computer  program  was  written  to  calculate  the 

.  1  •  ’  ■  i 

collision  probability  for  two  aircraft  making  the  simultaneous 

,  i 

curved  approaches  to  parallel  runways  shown  in  Fig.  5.  One 

;  i  i  ■  1 

aircraft  was  assumed. to  be  relatively  fast  (approach  speed 

>i  * 

\ 

120  knots),  the  other  relatively  slow  (approach  speed  80 
1  knots) .  The  details  of  the  approaches  are  given  in  tfabld  1 

and  Fig.  5.  The  aircraft  volumes  were  taken  to  be  right 

*  1  • 

circular  cylinders . 

!  ) 

!  The  statistics  were  assume^  to  be  Gaussian  and  the 
»  , 

standard  deviations  of  position  in  the  y  and  z  directions 

were  assumed  to  diminish  as  they  approached  the  runway  (see 

\ 

*  *  i 

Fig.  6).  This  'improved  accuracy  occurs,  in  practice,'  for 

1  .  .  ‘i  • 

r  ‘  1  ,  < 
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TABLE  1.-  DETAILS  OP  CURVED  APPROACHES 


Vehicle  1 

Vehicle  2 

Units 

Path  Segment  1 

duration 

315. 

60. 

sec 

ground  speed 

133. 

200. 

ft/sec 

descent  rate 

10. 

13.3 

ft/sec 

initial  altitude 

6600. 

9200. 

ft 

Path  Segment  2 

duration 

45. 

180. 

sec 

ground  speed 

133. 

200. 

ft/sec 

descent  rate 

10. 

13.3 

ft/sec 

initial  altitude 

3450. 

8400. 

ft 

Path  Segment  3 

duration 

450. 

600. 

sec 

ground  speed 

133. 

200. 

ft/sec 

descent  rate 

6.7 

10. 

ft/sec 

initial  altitude 

3000. 

6000. 

ft 

Path  Segment  4 

duration 

90. 

60. 

sec 

deceleration 

1,48 

3.33 

ft/sec/sec 

I  descent  rate 

0. 

0. 

ft/sec 

j  initial  altitude 

0. 

0. 

ft 
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VEHICLE  2 


120  KTS 


SEGMENT  I ' 


SEGMENT  2- 


- ► 


SEGMENT  4 


SEGMENT  3 


RUNWAYS(  0  ■  distance  between  runway  centerline* 


SEGMENT  4  SEGMENT  3 


■SEGMENT  2 


SEGMENT  I 


VEHICLE  I  Jl  80  KTS 


GEOMETRY  OF  VEHICLES 

VEHICLE  I :  RIGHT  CIRCULAR 
CYLINDER, HEIGHT  10  FT., 
DIA.50  FT. 

VEHICLE  2:  RIGHT  CIRCULAR 
CYLINDER, HEIGHT  30  FT. 

DIA.  150  FT. 


Figure  5.  Geometry  of  Simultaneous  Curved  Approaches, 
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NOTE  (1)  VEHICLE  #1-Y  POSITION-USE  SCALE  1 
VEHICLE  #2-Y  POSITION-USE  SCALE  1 
VEHICLE  #1-Z  POSITION-USE  SCALE  2 
VEHICLE  #2-Z  POSITION-USE  SCALE  3 

NOTE  (2)  THE  STANDARD  DEVIATIONS  OF  VELOCITY 
ERROR  ARE  CONSTANT  AT  12FT/SEC  IN 
THE  X  AND  Y  DIRECTIONS  AND  5.3  FT/SEC 
IN  THE  Z  DIRECTION 


Figure  6.  Standard  Deviations  of  Position  Vs. 
Distance  from  Touchdown. 


several  reasons:  The  ILS  errors  are  approximately  angular 
errors,  and  diminish  as  the  source  is  approached.  In  addition, 
aircraft  control  errors  diminish  with  length  of  time  on  the 
beam  and,  finally,  visual  cues  improve  as  touchdown  is 
approached.  None  of  these  circumstances,  however,  improves 
positional  guidance  accuracy  in  the  direction  of  the  runway, 
so  the  standard  deviations  of  the  errors  in  the  x  direction 
were  assumed  constant  at  6000  ft  for  each  vehicle. 

The  formula  for  the  large  relative  velocity  case  given 
by  (2)  was  found  to  be  adequately  accurate  even  though  the 
relative  velocity  passed  through  zero  shortly  after  the 
touchdowns.  This  accuracy  was  verified  by  comparing  the 
result  of  using  (2)  with  that  obtained  from  the  exact 
integration  of  (1) ,  which  does  not  contain  any  assumptions 
on  relative  velocity.  The  two  answers  agreed  to  within  8 
percent.  Since  the  computation  time  for  the  large  relative 
velocity  case  (2)  is  about  1/20  that  for  the  exact  integra¬ 
tion,  formula  (2)  was  used  throughout  the  runs.  It  should 
be  noted  that  formula  (2)  does  not  assume  independence  of 
position  and  velocity  statistics  and,  in  fact,  requires  no 
velocity  statistics  at  all.  Further,  although  Gaussian 
position  statistics  were  assumed,  any  other  statistics  may 
have  been  used  in  the  computer  runs  with  little  change  in 
the  program. 


RESULTS  OF  COMPUTER  RUNS 

The  results  of  the  computer  runs  are  shown  in  Pigs.  7 
through  9.  The  data  of  Figs.  5  and  6,  and  of  Table  1  apply 
to  all  runs  except  as  noted  in  the  following  discussions. 

Fig.  7;  (Variation  of  collision  probability  with  time 

between  touchdown.)  If  the  two  aircraft  come  abreast  at  the 

touchdown  point  (T«0)  the  collision  probability  is  very  small, 
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less  than  10  .  As  T  increases,  the  slower  aircraft  (vehicle 

1)  touches  down  later  than  the  faster  aircraft  (vehicle  2) 
and  the  point  of  passing  occurs  farther  up  the  glide  slope. 

Since  the  lateral  position  errors  increase  with  distance  from 
touchdown,  (Fig.  6) ,  the  collision  probability  increases  as 
the  point  of  passage  moves  away  from  the  touchdown  point. 

The  sharp  drop  in  collision  probability  beyond  T  =  100  is  due 
to  the  geometry  of  the  turn-on  for  vehicle  1. 

In  order  to  obtain  conservative  estimates,  the  runs  of 
Fig,  8  and  Fig.  9  were  made  with  T  =  80,  which  maximizes  the 
collision  probability.  If  the  runways  were  truly  independently 
operated,  the  times  between  touchdowns  would  be  distributed 
uniformly  between  0  and  the  time,  t,  between  successive  landings 
on  the  same  runway.  Since  t  often  is  about  120  seconds,  the 
assumption  is  not  overly  pessimistic. 

Fig.  8:  (Variation  of  collision  probability  with  runway 
separation.)  The  Figure  shows  that  cutting  the  runway  separa¬ 
tion  from  5000  ft  to  2500  ft  increases  the  probability  of  a 
collision  by  a  factor  of  6000,  other  things  being  equal. 
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PROBABILITY  OF  COLLISION 


Figure  7.  Probability  of  Collision  Vs.  Time 
Between  Touchdowns 
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This  result  is  reassuring  or  alarming,  depending  on  one’s 
view  of  the  safety  of  the  present  5000-ft  separation  and  one's 
belief  in  the  applicability  of  Gaussian  statistics  to  the 
problem. 

Fig.  9:  (Probability  of  collision  vs  lateral  navigation 
error.)  The  lateral,  or  y,  component  of  position  accuracy 
was  assumed  to  be  the  same  for  both  aircraft  and  was  chosen 
to  be  750  ft  for  the  runs  of  Figs.  7  and  8.  In  the  present 
case,  however,  it  is  allowed  to  vary.  Fig.  9  shows  that  the 
collision  risk  drops  drastically  when  lateral  placement  errors 
go  below  750  or  500  ft.  This  suggests  a  trade-off  between 
the  cost  of  lateral  guidance  equipment  and  the  cost  of  expanded 
airports . 
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CONCLUSION 

The  formulas  derived  allow  the  calculation  of  collision 
probability  between  two  aircraft  flying  arbitrary  curvilinear 
paths.  They  represent  a  generalization  of  the  statistical- 
probabilistic  method  of  analysis ,  which  heretofore  has  been 
available  only  for  constant  altitude,  rectilinear  paths.  The 
analysis  results  in  a  general  integral  expression  for  the 
collision  probability;  simpler  integral  formulas  are  derived 
for  the  cases  of  large  relative  velocity,  zero  relative  ve¬ 
locity,  and  a  spherical  collision  surface.  In  only  the  simplest 
Gaussian  cases  is  an  explicit  formula  available.  Because  high 
accuracy  is  not  essential  it  is  possible  to  compute  by  machine 
the  appropriate  integral  expression  formula  (2)  for  practical 
cases  such  as  curved  landing  and  takeoff  paths. 

The  computer  calculation  of  formula  (2)  for  curved 
approaches  to  parallel  runways  shows  how  the  collision  proba¬ 
bility  varies  with  time  between  touchdowns,  runway  separation, 
and  lateral  position  error. 
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